Scanning SQUID microscopy of vortex clusters in multiband superconductors 
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In type-1.5 superconductors, vortices emerge in clusters, which grow in size with increasing mag- 
netic field. These novel vortex clusters and their field dependence are directly visualized by scanning 
SQUID microscopy at very low vortex densities in MgE>2 single crystals. Our observations are elu- 
cidated by simulations based on a two-gap Ginzburg-Landau theory in the type-1.5 regime. 
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In magnesium diborideA, MgB2, the superconducting 
gaps open for both the two-dimensional (2D) cr-band 
(gap size A ff =7.1 meV- 2 -^) and the 3D 7r-band (A 7r =2.2 
meV— 4). The possibility of type-1.5 superconductivity 5 
has been suggested for clean single crystals of MgB 2 , 
which lie in the type-1.5 regime: K a =3.7> l/y/2 (type- 
2) and K 7r =0.66< l/v2 (type-1), where n a (a — a, if) is 
the Ginzburg-Landau (GL) parameter for each band esti- 
mated from the band structure calculations^. In contrast 
to conventional type-2 superconductors^, in type-1.5 su- 
perconductors, vortex stripes and gossamer-like vortex 
patterns, i.e., vortex clusters emerge at relatively low 
applied fields, 5 which is due to a competition between 
attractive (type-1) and repulsive (type-2) vortex inter- 
actions governed by a two-gap GL theory^ 12 . Interest- 
ingly, a substantial difference in vortex structure between 
type-1.5 and type-2 superconductors, as a fingerprint of 
type-1.5 superconductivity, is expected at very low vor- 
tex densities because the intervortex distances in clusters 
are likely to be almost independent of the applied field 
in type-1.5 superconductors, whereas the intervortex dis- 
tances in type-2 superconductors follow the conventional 
dependence (</> /-B) 1/2 , where o (=2.O7x KT 15 T m) is 
the flux quantum and B the magnetic field. This mo- 
tivates strongly direct vortex visualization experiments 
aimed at investigating vortex structure at very low vor- 
tex densities in MgB2 single crystals. Additionally, it 
is important to verify the existence of type-1.5 supercon- 
ductivity in high quality MgB2 crystals different from the 
ones used in Ref. [a. 

In order to study vortex structure in MgB2, we made 
scans with a scanning superconducting quantum inter- 
ference device (SQUID) microscope on single crystals 
of MgB2 which were grown by a pressure synthesis 
technique^ (a crystal grown by this technique was used 
elsewhere^). The superconducting transition temper- 
atures T c of the crystals measured by a SQUID mag- 
netometer are 38.5 K [the transition width AT C =0.8 K 
(10-90 % criterion)]. The surface of crystals was con- 
firmed not to have any cracks or holes by using a field 
emission scanning electron microscope in the secondary 
electron regime. High-resolution transmission electron 



microscope images and electron diffraction patterns for 
the crystals show that the crystals have no grain bound- 
aries. No impurity contamination was detected within 
0.1% accuracy by an electron probe microanalyzer. 

Our scanning SQUID microscope is based on prototype 
SQM2000±S, (SII NanoTechnology Inc.). We fabricated a 
SQUID sensor, which is a SQUID magnetometer linked 
with a circular Nb pickup loop. The Nb pickup loop with 
an inner diameter of 8 /im scans the surface of a sample, 
keeping a few micrometers away from the surface. A 
magnetic flux sensitivity is 5 /i^o/Hz 1 / 2 . A field noise 
referred to the pickup loop is ^10~ 10 T/Hz 1 / 2 . A spatial 
resolution is ^2.5 /im. Before scans, as-grown crystals 
were cooled from a temperature higher than T c to 4.2 K 
under the magnetic field parallel to the c-axis and then 
scans were made on the afc-plane at 4.2 K. Scanning steps 
are 0.5 //m and 1.0 /im for MgB2 single crystals and a 
Nb film, respectively. 

We model vortex systems in a type-1.5 superconduc- 
tor in the same way as in Ref. [5j. We invoke the fol- 
lowing two-gap GL theory to derive the vortex-vortex 
interaction numerically. A two-gap GL free energy 
functional^— is the sum of two single-band GL func- 
tional with the Josephson coupling term corresponding 
to the interaction between two bands: 



F = 



J dr 3 (F a 



F n + F„ n + (1) 
2^o 



where F a (a=a, n) is the free energy of each band and 
h (=|h|) is a magnetic field (h= rot A, where A is the 
vector potential) 

27T 

F a = 2E a \^ Q \ 2 + \E a \\^ a \ i + C\{-tV + — A)* Q | 2 , (2) 



where E a is the 
C=(0 o /27r) v /|£; a |/ At oK3 
term: F n7r 



condensation energy- and 
F ar: is the Josephson coupling 



, , ,. , T . , , where E 1 is the 

coupling energy, where we take a value used in Ref. 
ITtI The vortex-vortex interaction energy between 
two vortices is the sum of increased energies of each 
vortex caused by the presence of the other. It can be 



2 



Magnetic field (^T) 
10 20 30 40 50 



1 20 30 40 50 



< 3 k • 




Position ((jm) 
5 10 15 20 25 



Magnetic field (|iT) 

10 20 30 40 50 10 20 30 40 50 50 



100 





(c > o Data t 




10 20 30 
Distance (|im) 



FIG. 1: (color online) Scanning SQUID microscope images 
of vortices at (a) 1 /iT and (b) 10 /xT for a MgB2 single 
crystal, the dimensions of which are approximately 80x60 
/im 2 . Dotted lines show the edges of the sample. A scale 
bar corresponds to 10 /im. The integral of the magnetic field 
within a solid line in (a) gives the flux ~<^o- (c) The cross 
section of the vortex image along a broken line in (a) (open 
circles) and the model fit (solid line) (see text), (d) Scanning 
SQUID microscope image of vortices at 10 fiT in a Nb film 
with a thickness of 200 nm (T c =9.1 K). (e) The 2D multi- 
fits of the SQUID image in (b) to numerical calculations by a 
London model (see text), (f) and (g) The distribution of the 
intervortex distance derived from (e) and (d), respectively. 
A solid line in (g) represents a fit to the Gaussian function. 
An arrow indicates the intervortex distance calculated from 
Wo/B) 1/2 . 



derived by numerically minimizing the two-gap GL free 
energy of two vortices with a variational procedure 18 . 
The minimization gives the interaction consisting of 
short-range repulsion and weak long-range attraction, 
which is similar to that reported in Ref. |l8|. 

An equation of motion for a vortex i is Fj = F™ + 
Fj = ijVi, where F™ is the vortex- vortex interaction 
term we calculate above, Fj the thermal stochastic force 
which satisfies (Ff (i))=0 and (Ff (t)Fj(f)) = 2r]8 ij 5{t- 
t')fcsT, and rj the viscosity, where we take a value esti- 
mated in Ref. [B]. The systems with no pinning cen- 
ters were initially prepared in a high temperature molten 
state and then annealed with two million temperature 
stepsJ^ We made the system stable during 2000 time 
steps in each step of temperature. 

Figures 1(a) and 1(b) show a scanning SQUID micro- 
scope images of vortices at very low vortex densities for 
a MgB2 single crystal. The size of a vortex in the image 
does not correspond straightforwardly to the penetration 
depth A but rather to a combination of A and the extent 
of a stray field emanating from a vortex at the sample 
surface. In order to estimate A it is necessary to compare 
the cross section of the vortex image with appropriate 
numerical calculations where the extent of the stray field 
is taken into account. 15 We calculate numerically the spa- 
tial distribution of a magnetic field from a vortex above 
the sample surface, using a London mode l 15 ' 20 ' 21 . 
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FIG. 2: (color online) Scanning SQUID microscope images of 
vortices at (a) 10 itT, (b) 20 f/F, and (c) 50 itT for a MgB 2 
single crystal. A scale bar corresponds to 10 /im. (d)-(f) The 
2D multi-fits of the SQUID images in (a)-(c) by a London 
model (see text), (g)-(i) The locations of vortices in (d)-(f) 
obtained by the 2D multi-fits of the images. Broken lines 
show the edges of the sample. 



Figure 1(c) shows a fit of the cross section of a vortex 
image in Fig. 1(a) by the numerical simulations. We 
obtain A=0.01-0.13 /im (z=1.7±0.1 /im, where z is the 
distance between the sample surface and the pickup loop) 
from the fit. The z value was determined by the fit of a 
vortex image in a Nb film (known parameter A=50 nm) 
mounted on a sample holder with a MgB2 crystal. The 
range of A is a residual based on the error bar in the 
z value which was assigned by doubling of the \ value 
in the fit in a Nb film (a doubling of the x value only 
gives statistical errors). Due to the error it is difficult 
to specify how close our A value is to the band structure 
calculations 6 . However, first-principles calculations^ 2 ^ 
of the penetration depth in the clean limit fall on our A 
value within the error bar. 

In type-2 superconductors with the weak pinning, the 
Abrikosov lattice is often seen in field cooling experi- 
ments even if a vortex density is quite lowJ^ the vortex 
lattice made of vortices strongly interacting with each 
other at temperatures very close to T c is frozen with de- 
creasing temperature. Vortex patterns at very low vortex 
densities in a MgB2 single crystal, which has the weak 
pinningfil seem to be different from that in the conven- 
tional weak pinning system. At B—l /xT vortices are 
located far away from each other. [Fig. 1(a)] However, 
with increasing applied field vortices come closer to each 
other to form compact groups as shown in Figs. 1(b) 
and 1(e), which is never seen in superconductors with 
the weak pinning, such as MoGe films^ and NbSe2 sin- 
gle crystals 5 -. 
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FIG. 3: Scanning SQUID microscope images of vortices at 20 
fiT obtained from different cool downs. The left image is the 
same as Fig. 2(b). 



In a Nb film with the strong pinning, disordered vortex 
patterns are observed at low vortex densities, as shown 
in Fig. 1(d). We now compare qualitatively these vortex 
patterns. Figure 1(e) shows the 2D multi-fits^ of the 
scanning SQUID image in Fig. 1(b) to numerical cal- 
culations from a London model2&2i. Fits were made by 
parameterizing the coordinates of each position of vor- 
tices with z and A determined above. From the fits it 
turns out that the number of vortices is 23, which is in- 
deed consistent with 23 estimated by SB/(f> , where S is 
the areas of a sample. The distribution of the intervortex 
distance can be derived from the fits, as shown in Fig. 
1(f). This distribution is a bimodal and it clearly splits 
into two: one ranging from ^2 /jm to ^7 (im and another 
from ~15 /im to ^30 /jm. In principle, the former corre- 
sponds to the intragroup distribution and the latter cor- 
responds to that concerning the intergroup. On the other 
hand, the vortex distribution in a Nb film can be fitted by 
the Gaussian function, as shown in Fig. 1(g). The cen- 
ter of the peak in the Gaussian function is at 14.5 /jja, 
indicating that the average intervortex distance is 14.5 
fxm. This agrees with the intervortex distance (fio/B) 1 / 2 
of ^14 /jm estimated by assuming the presence of the 
Abrikosov lattice. This fact indicates that the vortex 
pattern in Fig. 1(b) is not a disordered Abrikosov lattice 
but vortex cluster array. On the other hand, the vor- 
tex pattern in Fig. 1(d) does correspond to the strongly 
disordered Abrikosov lattice. 

Figures 2(a)-2(c) show the variation of the vortex 
structure with applied fields. After vortices form local 
clusters in Fig. 2(a), some islands of magnetic flux and 
the semi-Meissner regions come out of the grouping with 
increasing applied field. [Figs. 2(b) and 2(c)] The posi- 
tions of the vortices can be determined by the 2D multi- 
fits [Figs. 2(d)-2(f)] of Figs 2(a)-2(c), as shown in Figs. 
2(g)-2(i). The position maps show that vortices line up 
locally in the early stage [Fig. 2(g)] and then small vor- 
tex clusters merge with each other, [Fig. 2(h)] to grow 
into bigger clusters with applied fields. [Fig. 2(i)] The 
number of vortices changes to 49 in Fig. 2(h) and then 
114 in Fig. 2(i). These are reasonably well approximated 
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FIG. 4: Simulations of systems consisting of (a) 6, (b) 27, 
(c) 122, and (d) 500 vortices which have the interaction en- 
ergy calculated from a two-gap GL theory for the type-1.5 
condition (see text). A scale bar corresponds to 10 /im. 



by 46 and 116 estimated from SB /fa, respectively. 

We note that the vortex array varies considerably 
in the locations of clusters measured by different cool 
downs, while keeping the same general features, as seen 
in Fig 3. Since the pinning sites are fixed, this implies 
that the vortex clustering is not due to the presence of 
pinning centers but is rather to be attributed to the pe- 
culiar vortex- vortex interaction combining repulsion and 
attraction. 

Figures 4(a)-4(d) show results of the simulation men- 
tioned above. In type-1.5 superconductors vortices begin 
to form clusters consisting of a few vortices at very low 
vortex densities, [Fig. 4(a)] which does not coincide with 
the experiment displayed in Fig. 1(a). Perhaps, the exis- 
tence of some pinning centers may be responsible for this 
fact. This will be considered in the future work. The 
number of small clusters increases with applied fields, 
[Fig. 4(b)] which apparently corresponds to Fig. 2(g). 
Then small clusters start to merge with each other, [Fig. 
4(c)] to grow into bigger clusters with a further increase 
of an applied field [Fig. 4(d)]. This simulation result co- 
incides quite well with experimental data [Figs. 2(h) and 
2(i)]- 

The simulations show that in a vortex cluster the dis- 
tances between vortices and their nearest neighbors are 
almost independent of the applied field (~1 fim). There- 
fore, these distances are not determined by the expres- 
sion (0o / B) 1 / 2 but rather by the competition between 
short-range repulsive and long-range attractive interac- 
tions between vortices. In Fig. 5 the average intervortex 
distance in clusters obtained from scanning SQUID and 
decoration images^ is shown as a function of applied field, 
compared with the dependence {ifio/B) 1 / 2 . The average 
distances indeed do not follow any dependence close to 
(4>q/ B) 1 / 2 but follow instead a very weak linear field de- 
pendence. The intervortex distance in clusters, indepen- 
dent of the applied field, can be identified as a unique 
property of a system of vortices with short-range repul- 
sion and long-range attraction, i.e., vortices induced in 
a type-1.5 superconductor. This fact indicates that vor- 
tex clusters definitely emerge even in crystals grown by 
a different crystal maker, implying that in general vor- 
tices tend to form vortex clusters in high quality single 
crystals of MgB 2 £ 

A difference between (</>o/-B) 1//2 and the intervortex 
distance in clusters decreases as an applied field increase 



4 




o b i i i i i t 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 
Applied field (mT) 



FIG. 5: The average intervortex distance in vortex clusters 
as a function of applied field. These distances are obtained 
from fits of scanning SQUID images in Figs. 2(g)-2(i) (filled 
circles) and decoration images* at 0.1 mT and 0.5 mT (filled 
squares). Error bars indicate the range of the distribution 
of vortex distances in clusters. A solid line represents the 
dependence (^o/-B) 1//2 . 

and these two length scales become equivalent around 0.6 
mT. Any vortex cluster is unlikely to be observed above 
0.6 mT. In fact Moshchalkov et al£ have reported no 
clearly visible vortex clusters at 1.0 mT. The observations 



of long-range Abrikosov lattices at 20 mT and 0.2 T— ^ 
are consistent with this scenario. Although Vinnikov et 
ati^i have obtained a decoration image of a disordered 
vortex lattice at 0.44 mT, it would be too difficult to 
find vortex clusters, a size of which is more than 20 /im 
around 0.5 mT, as reported in Ref. 

In conclusion, scanning SQUID microscopy shows that 
in high quality MgB2 single crystals vortices line up lo- 
cally at very low vortex densities and grow into vortex 
clusters with increasing applied field. The average dis- 
tances between nearest neighboring vortices forming clus- 
ters are practically independent of the applied field B, in 
contrast to the conventional (1/ B) 1 / 2 behavior. The sim- 
ulations of a vortex system of a type- 1.5 superconductor 
indeed reproduce these observations. Our experimental 
data, as well as consistency between experiment and the 
simulations, clearly demonstrates that high quality MgB2 
single crystals lie in the type-1.5 regime. 
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